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We study here different regions in phase diagrams of the spin- 1/2, spin-1 and spin-3/2 one dimen- 
sional antiferromagnetic Heisenberg systems with frustration (next-nearest-neighbor interaction J2) 
and dimerization (S). In particular, we analyze the behaviors of the bipartite entanglement entropy 
and fidelity at the gapless to gapped phase transitions and across the lines separating different phases 
in the J2 — 8 plane. All the calculations in this work are based on numerical exact diagonalizations 
of finite systems. 
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I. INTRODUCTION 



Matter can appear in different quantum phases with 
exotic properties like charge density wave, magnetism, 
superconductivity, and so on. Studies of these phases 
and the transitions from one phase to the other are im- 
portant and interesting for both academic and techno- 
logical reasons. Two major tools from quantum infor- 
mation theory have been used extensively in recent years 
for studying the quantum phases of a system: quantum 
entanglement and fidelity. The idea of quantum entangle- 
ment originated in the study of quantum correlations of 
many-body systems^. It is expected that even for mod- 
erately large system sizes, the entanglement entropy can 
identify the values of the parameters of the Hamiltonian 
where a quantum phase transition (QPT) occurs because 
the quantum correlations of the systems change signifi- 
cantly occurs when one goes across such a transition 3 -. 
In recent years, entanglement entropy has been used to 
study quantum critical regions in various systems^—. 
Quantum fidelity is a measure of how little the ground 
state of a system changes as one changes the parameters 
of the Hamiltonian. A large change in the fidelity is an- 
ticipated close to a QPT even if the system size is not 
very largei^r— . 

In this paper, we study different quantum phases 
and quantum critical regions with frustration J2 (next- 
nearest-neighbor coupling) and dimerization S (an alter- 
nation in the nearest-neighbor couplings) of the spin-1/2, 
1 and 3/2 Heisenberg antiferromagnetic systems in one 



dimension by calculating the entanglement entropy and 
fidelity in the ground state of the system. In particu- 
lar, we study the transition from a gapless to a gapped 
phase and the changes in the spin structure across differ- 
ent phase lines in the J 2 — 5 plane for the spin-1/2, 1 and 
3/2 systems. This numerical study helps us to estimate 
the locations of different critical points and the lines sep- 
arating different phases occurring in the phase diagrams 
of these systems. We also compare our results with those 
reported previously (using different techniques) whenever 
possible. 

There are some related works which we briefly men- 
tion here. The bipartite entanglement for the spin-1/2 
J\ — J2 model (without dimerization) has been studied 
in different contexts, like the transition from the Neel to 
the spiral phase and the gapless to gapped phase transi- 
tion along the J2-axis (see, for example ; 29 ' 30 ). The gap- 
less to gapped phase transition has also been studied by 
the fidelity of the first excited state (as a function of J2) 
with periodic boundary conditions^!. The role of entan- 
glement between distant sites has been studied for this 
model with frustration and dimerization^ 2 -. To the best 
of our knowledge, the phase diagram of the spin-1 and 
spin-3/2 J\ — J2 model (with or without dimerization) 
has not yet been studied using entanglement entropy and 
fidelity. In this paper we report results on the entangle- 
ment entropy and fidelity of the J\ — J2 — 5 model over 
the J2 — 5 plane. 

Our paper is organized in the following way. In section 
2, we discuss how to calculate the entanglement entropy 
and fidelity. We then introduce the Hamiltonian that we 
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FIG. 1. Phase diagrams of the (a) spin-1/2 and (b) spin-1 chains in the J2 — 8 plane. 



study in section 3. In section 4, we give a brief intro- 
duction to the numerical techniques we use in our work. 
We then present the entropy and fidelity results for the 
spin-1/2, spin-1 and spin-3/2 systems in section 5. We 
conclude our paper in section 6. 



II. ENTANGLEMENT ENTROPY AND 
FIDELITY 

A pure state of a bipartite entangled system can be 
written as \tp) = ^^Cij\<j)i) l \cj>j) r , where | and \(pj) r 

ij 

are the basis states of the left and right blocks respec- 
tively. The reduced density matrix (RDM) of the left 
block, pi = Tr r (\ip){ip\), is calculated by tracing out the 
degrees of freedom of the right block. The elements of 
the RDM pi are given by 

Pij = ^CikCjk. (1) 

k 

The von Neumann entropy of a block is given by S = 
—Tr{p logo, p) or 

S = -J2 \ log 2 X u (2) 

i 

where the Aj's arc the eigenvalues of p. 

Fidelity measures how little a particular wave function 
(for instance, the ground state) changes with the param- 
eters of a model Hamiltonian. This is quantified by the 
overlap of the wave function at two different parameter 
values. If p is a parameter then the fidelity is given by 

F = |(*(p)|tt(p + a)>|, (3) 



where a is a small variation in p. In our case, both J 2 
and 6 are parameters with respect to which we have cal- 
culated fidelity and we have taken the change in the pa- 
rameter to be 10 -2 in the numerical calculations. 



III. DESCRIPTION OF THE SPIN MODEL 

We study the Heisenberg Hamiltonian for the antifcr- 
romagnctic chain with both nearest-neighbor and next- 
nearest-neighbor couplings and dimerizatio n 33 ' 34 . We 
will use this Hamiltonian for the spin-1/2, spin-1 and 
spin-3/2 systems; it is given by 

2JV-1 

h = j 1 y / (i-(-ir s)Si -s i+1 
i=i 

2N-2 

+J 2 ]T S t ■ S i+2 , (4) 
1=1 

where J\ is the nearest-neighbor interaction (we take 
Ji=l for our study), S (0 < S < 1) is the dimerization and 
J2 (0 < J2 < 2) is the next-nearest-ncighbor interaction. 
In our entropy contour plot we have taken the range of 
j r 2 from to 2. 

The phase diagrams of the spin-1/2 and spin-1 chains 
are different from each other in the J2 — 5 plane (see 
Fig. Q]). The spin-1/2 system undergoes a QPT from a 
gapless phase to a gapped phase at J 2c = 0.2411 ±0.0001 
without dimerization (S = 0)2£, while the rest of the 
phase diagram is gapped. The line 2 J 2 + 5=1 separates 
the Neel phase (region A) from the spiral phase (region 
Bp.. 

There are a number of distinct phases in the spin-1 
system in the J2 — 5 plane. Region I represents the spin 
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FIG. 2. Some legal valence bond diagrams of a spin-1/2 chain 
with 8 sites and total spin S = 0. 



Peierls gapped phase. Regions II and III are Haldanc 
gapped and spiral gapless regions respectively. In these 
two regions the ground state of the open chain is four- 
fold degenerate. Region IV is a spiral gapped phase. A 
gapless phase exists along the critical line 'a' that lies 
between (0, 0.25) to (0.22±0.02, 0.20±0.02) in the J 2 - 5 
pian o 36 ' 37 . A line 'c' separating the regions II and III 
extends from (0.73, 0) to (0.65, 0.05); on that line the gap 
appears to be zero (to numerical accuracy). Along the 
dotted line 'b' the gap is finite. The line 2 J 2 +<5 = 1 starts 
at the point P = (0.432, 0.136) and extends up to (0, 1), 
separating regions I and IV. Another line ('d') starts at 
(0.39, 0) and ends at point 'P', separating regions II and 
III. 

The phase diagram of spin-3/2 system has not been 
studied yet. We will show below that the entanglement 
entropy and fidelity of this system can give some insights 
into its phase diagram. 



IV. NUMERICAL TECHNIQUES 

For our calculations, we use the M s basis (eigenstates 
of the z component of total spin) 3 ^. These basis states 
are orthonormal and it is easy to obtain the RDM (which 
is used to calculate the entanglement entropy) when the 
states of the system are expressed in this basis. On 
the other hand, most of the results can be understood 
qualitatively using valence bond (VB) theory 3 ^. In this 
theor y 40 ' 41 , the basis states in the singlet space are ex- 
pressed as products of pairwise singlets, which follow the 
Rumer-Pauling rules, to avoid overcomplctcncss of the 
VB states. Fig. shows some VB-diagrams of a spin- 
1/2 chain with 8 sites and total spin 5 = 0. The VB 
state (a) is a Kekule state which is a product of nearest- 
neighbor singlets. 

Since our Hamiltonian conserves total spin, all its 
eigenstates are also eigenstates of total spin. There- 
fore, the eigenstates that we obtain by diagonalizing the 
Hamiltonian in the constant M s basis will be linear com- 
binations of the VB basis states. It may be worth men- 
tioning here that appropriate linear combinations of con- 
stant M s basis states can give different VB basis states^. 
Now, a VB basis state contributes to the bipartite entan- 
glement entropy of a state under study if the boundary 
between the two blocks of the system cuts a singlet line. 
If the boundary does not cut a singlet line, its contribu- 
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FIG. 3. Ground state entanglement entropy of the spin-1/2 
system with different chain lengths at S — 0. The lower set 
of plots for 4n systems and the upper set of plots is for 4n+2 
systems. 



tion is assumed to be zero. For example, if the boundary 
goes through the sites 4 and 5 of the system, the en- 
tropy contribution of diagram (a) will be zero while that 
of diagram (b) will be be non-zero (Fig. [2]). Depending 
upon the entropy contributions of the VB basis states 
and their relative weights in a state under study, we can 
qualitatively understood the entropy of the stated. In 
generating the contour plot of entropy, we have calcu- 
lated the entropy over a grid of 201 J 2 values (0<J 2 <2) 
and 101 S values (0< S <1). 



V. RESULTS AND DISCUSSION 

In this section, we present numerical results for the 
entanglement entropy of finite size chains with two equal 
block sizes for the spin-1/2, 1 and 3/2. We also present 
results for the ground state fidelity of both the systems 
in the J 2 — 5 plane. 



A. The spin-1/2 system in J2 — S plane 

1. Uniform chain (8 = 0) 

The uniform spin-1/2 chain without dimerization (5 = 
0) goes through a gapless to gapped phase transition at 
J 2c — 0.2411 in the thermodynamic limit (N — > oo)2£, 
and its Neel phase is separated from the spiral phase at 
J 2 = 0.5. To study the behavior of the system around 
those points, we calculate the ground state fidelity and 
entanglement entropy of the system. 

We first consider the ground state entanglement en- 
tropy of finite size spin-1/2 chains with different J 2 val- 
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FIG. 4. Lowest and first excited state energies of 16 and 20 sites for different Ji values with total spin S=0 and 1 are shown in 
(a) and (b). Entropy of the second excited state (dashed line) and third excited state (solid line) of 16 and 20 sites for different 
Ji values in the M s = sector are shown in (c) and (d). 



ues. The bipartite entropy for different chain lengths 
(and equal block size) can be seen in Fig. [3] At Ji = 0.5, 
the entropy reaches a minimum; away from that point 
the entropy increases. For the systems with even block 
sizes, the minimum of the entropy goes to zero while for 
odd block sizes, this minimum is one. This result can be 
explained by noting that, at J 2 = 0.5 (the Majumdar- 
Ghosh pointed) , the ground state has Kekule state struc- 
ture (as in Fig. [210 • Depending upon the block size being 
odd or even, entropy will be finite or zero respectively 3 ^. 
As J 2 moves away from that point, the presence of other 
VB basis states (as in Fig. HJb)) in the ground state will 
become significant. Since these VB basis states will have 
a finite entropy contribution, the entropy of the ground 
state will increase as J2 moves away from the point. 

We do not observe any behavioral change in the ground 
state entropy at about J 2c- The first excited states in 
the singlet and triplet sectors cross as a function of J2. 
We have calculated the entropy of these two states as 
a function of J2. We find that the entropy of a state 
depends only on its spin and not on its energy. There- 
fore as a function of energy level ordering the entropy 
shows a jump at a value of J% which depends on the 
chain length. The jump in the value of the entropy can 
be seen from Fig. [4] In Figs. [4] (a) and (b), we see that 
the ground state (lowest singlet) and lowest triplet state 
arc non-degenerate for finite N; they become degenerate 
in the thermodynamic limi t 34 i 35 . However, the first ex- 
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FIG. 5. Convergence of the crossing points of the excited 
states for a spin- 1/2 chain. 

cited singlet and the first excited triplet states become 
degenerate near J2 C even for small values of N. In Fig. 
[5] we plot the J 2 (N) as a function of 1/ iV 2 and we see 
that J| extrapolates to J 2c in the thermodynamic limit. 
The extrapolated value of J| (0.2414) is very close to the 
reported value of J 2c = 0.24112^. 

Upon calculating the ground state fidelity (taking J 2 to 
be the variable parameter), we do not observe any bchav- 
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FIG. 6. Ground state entanglement entropy of a spin-1/2 chain with lengths of 6, 8, 10, 12, 14, 16, 18 and 20, for different 
values of 5 values; the line types for different lengths are the same as in Fig. [3] The upper set of curves is for chain length 
4n + 2 while the lower one is for chain length An (n being a positive integer). Near the minimum, the finite size dependence is 
very weak. 




FIG. 7. Ground state entanglement entropy contour plot of spin-1/2 chain with 20 sites, for J2 going from to 2 and S going 
from to 1. The entropy along the dashed line is very nearly zero and corresponds to a Kekule ground state. 



6 




ioral change at J 2 = 0.5 or J 2c . The fidelity of first excited 
states in both the singlet and the triplet sectors falls to 
zero near J2 C . The same thing has already been observed 
for the system with periodic boundary conditions^. 



2. Dvmerized chain (0 < 5 < 1) 

In the phase diagram of the spin- 1/2 system with 
dimerization, the Neel phase is separated from the spiral 
phase by the line 2J2 + 5 = 1. We study the ground 
state entanglement entropy of finite size chains in these 
phases with different S values (see Fig. |5J). We observe 
that the entropy of the system is minimum for values of 
J2 and 6 which fall on this line. For systems with even 
sized blocks, this minimum value is zero while for odd 
sized blocks the minimum is one. The reason for this is 
similar to the case S = as given earlier. 

We have calculated the ground state fidelity (with J2 as 
the variable parameter) of finite size systems in the J2 — 5 
plane; we do not observe any sudden change in fidelity 
along the 2 J2 + 5 = 1 line. This can be explained by the 
fact that the phases of the system on both sides of the 
line are gapped which implies that the ground state does 
not cross any excited state (i.e., the ground state does 
not change its character of being a singlet or a triplet) 
when we cross this line in the parameter space. 



3. Spin-1/2 entropy phase diagram (contour plot) 

We study the gapless to gapped phase transition and 
change in the spin structure (order-disorder change) 
along the 2J2 + 5 = 1 line using the entropy contour 
plot for a spin-1/2 chain with 20 sites (see Fig. [7|). The 
entropy is zero along this line (shown by dashed line in 
Fig. 0. The density of the entropy contour lines shows 
whether the system is in a gapless or gapped phase. From 
this figure we see that the gapless region between J2 = 
and J ic has a higher density of entropy contour lines. In 
the rest of the figure the density of the contour lines is 
lower which shows that the rest of the phase diagram is 
gapped. At higher values of J2, the spin-1/2 chain be- 
haves like two decoupled chains with a weak coupling ( J\) 
between them. The stronger interaction (J2) is respon- 
sible for the higher entropy. For higher values of both 
J2 and 5 values, the spin-1/2 chain behaves like a spin 
ladder with a gapped phase. 



B. The spin-1 system in J2 — S plane 

As mentioned in Sec. 3, there are many phases in the 
J2 — S phase diagram of the spin-1 system. To study these 
phases, we first calculate the bipartite entanglement en- 
tropy in the ground state for finite systems. We calculate 
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the ground state entropy for different values of 5 (with 
large 5) and for different chain lengths (see Fig. |8|). The 
entropy is zero along the line 2 J2 + 5 = 1 for the systems 
with even sized blocks and non-zero for systems with odd 
sized blocks just as in the spin-1/2 case. As in the spin- 
1/2 case, the finite size effects are very weak near the 
minima of the entropy (Fig. [5]). 



1. Spin-1 entropy phase diagram (contour plot) 

In the Fig. [IJb), we see that the ground state is four- 
fold degenerate in regions II and III for an open chain. 
Hence we used the spin parity symmetry to break the 
degeneracy between the states corresponding to the total 
spin S = and 1. Then we calculate the entanglement 
entropy for the lowest state in the even parity subspace 
(which is a singlet). We study the quantum phases and 
QPTs of the spin-1 chain with 16 sites using the contour 
plot of the ground state entropy in the J2 — 5 plane, 
where J2 goes from to 2 and 5 goes from to 1 (see 
Fig. |9]). We see in the figure that the line 2J2 + 5 = 1 
starts approximately at the point P = (0.432, 0.136) and 
extends up to (0, 1). Along this line the entropy is zero. 
Along the gapless lines 'a' and 'c' (in Fig. QJb)) the 
density of the entropy contour lines is higher, while the 
density is lower in the rest of the phase diagram. About 
the line 'd' (in Fig. Hfb)), we observe that the density of 




FIG. 10. Curve (i) shows the (J2,5 ca i) values (squares) cor- 
responding to the minimum entropy. Curve (ii) shows the 
(J2,5 ca i) values (circles) corresponding to the maximum ab- 
solute values of the first order derivative of the entropy w.r.t. 
S. Inset shows the points corresponding to minimum fidelity 
(taking 5 as the variable parameter). All the results are ob- 
tained for the 16 site spin-1 chain. 



contour lines is much lower compared to the regions near 
by 

In curve (i) of Fig. [TU1 we plot S ca i versus J2 for the 
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FIG. 11. Ground state entanglement entropy contour plot of a spin-3/2 chain with 12 sites, with J2 going from to 2 and 5 
going from to 1. Dashed line corresponds to the contour along which the entropy is very nearly zero. 



points corresponding to the minimum value of the en- 
tropy in the J2 — S plane; this curve goes from the point 
'P' to the point (0, 1). This curve is seen to follow the line 
2J-2 + S = 1 as expected; we see a step staircase instead of 
a straight line because we have calculated the entropy for 
discrete values of J 2 and 5. To detect the gaplcss phase 
along the lines 'a' and 'c', we calculate the first order 
derivative of the entropy (using the three-point differen- 
tiation formula) along the <5-axis for different values of J2 
in this plane. In curve (ii) of Fig. 1101 we plot S ca i versus 
J2 corresponding to the maximum of the absolute value 
of the derivative. This curve follows the lines 'a' and 'c' 
closely. 



2. Spin-1 chain fidelity 

We calculate the ground state fidelity of spin-1 chains 
with finite sizes in the J2 — S plane. For small values of 
5, the ground states of finite size systems have multiple 
energy level crossings with the excited states. Because of 
these finite size effects, the fidelity of the ground states is 
not a reliable tool for studying phases in regions II and 
III. However in the regions I and IV, we find no energy 
level crossings in the ground states of finite systems. We 
calculate the ground state fidelity for systems with 6, 8, 
10, 12 and 16 sites for different values of S along the 
J2 axis in this plane. We find no sudden changes in the 
fidelity We show the plot of J2 versus S ca i corresponding 



to the minimum fidelity for a chain of 16 sites (see inset 
of Fig. [T0|) . This curve follows the lines 'a', 'b' and 'c' 
qualitatively separating regions II and III from regions I 
and IV. 



C. The spin-3/2 system in J2 — S plane 

In this section we study the phase diagram of the spin- 
3/2 system in the J 2 — <5 plane. To the best of our knowl- 
edge, the quantum phases of the Heisenberg spin-3/2 an- 
tiferromagnetic chain in the J2 — S plane has not been 
studied earlier. But from the field theory analysis of the 
spin chain, half-odd integer systems are gapless at S = 
and for small values of J2. With dimcrization (8^0), 
it is predicted that spin-3/2 system should be gapless at 
£ = 2/3 for J 2 = 02L4&. 



1. Spin-3/2 entropy phase diagram (contour plot) 

Wc study different quantum phases of spin-3/2 chain. 
We use spin parity symmetry to break the degeneracy 
(within numerical accuracy) of the ground state of this 
system. For 12 site chain, a contour plot of the entropy 
is shown in Fig. [TT] As in the spin- 1/2 and spin-1 cases, 
wc observe an order-disorder transition along 2J2 + 5 = 1 
line. The line starts approximately from the point Q = 
(J 2 = 0.38, 5 = 0.24) and extends up to (J 2 = 0, S = 1). 
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J 2 

FIG. 12. Curve (i) shows the (J2,S ca i) values (squares) cor- 
responding to the minimum entropy. Curve (ii) shows the 
(J2, 5 C ai) values corresponding to the maximum absolute val- 
ues of the first order derivative of the entropy w.r.t. S. Inset 
shows the points corresponding to minimum fidelity (taking 
S as the variable parameter). All the results are obtained for 
the 12 site spin-3/2 chain. 



For small J 2 values near S = in the contour plot, the 
line pattern is similar to that for spin- 1/2 case and quite 
different from the spin-1 case. This suggests that there 
can be a gapless phase at 5 — as predicted. For J 2 = 0, 
the line density is very high between 6 = 0.4 and 0.5; this 
suggests another gapless phase in this region as predicted 
by field theory. As in the spin-1 case, the density of lines 
is high at larger J2 values (about J2 = 1). This suggests 
a numerically gapless phase in this region. We also con- 
firm these gapless phases by comparing numerical energy 
gaps in those regions with that of a phase at large values 
of J2 and S where the density of lines is very low. For 
better understanding of these quantum phases we plot 
Seal versus J 2 values corresponding to the minimum en- 
tropy (curve (i) from Fig. [TO]) above the 'Q' point. This 
curve follows the 2 J 2 + S = 1 line similar to the spin-1/2 
and spin-1 cases. We also plot the points corresponding 
to the maximum absolute values of the first order deriva- 
tive of entropy with respect to 5; this is shown as curve 
(ii) in Fig. 1121 This is similar to the curve representing a 



numerically gapless phase for spin-1 system (curve (ii) in 
Fig. [TO]). This suggests that there can also be a gapless 
region along the curve (ii) for spin-3/2 system. The gap- 
less point ( curve (ii) of Fig. [TO]) at 6 = 0.45 at J 2 = is 
consistent with the value 5 = 0.431 reported^ 4 -. Note that 
this value is different from the field theory prediction of 

We calculate ground state fidelity of the 12 site spin- 
3/2 chain along different J 2 values with 5 as the variable 
parameter. We plot J 2 versus 5 correspond to minimum 
fidelity in this plane (see inset of Fig. [TO]) . The curve 
approximately follows the curve (ii) in Fig. 1121 
VI. CONCLUSION 

We have used entanglement entropy and fidelity as 
tools to study the different quantum phases and quan- 
tum critical regions of the spin-1/2, 1 and 3/2 chains in 
the J2 — S plane. For this study, we have employed ex- 
tensive exact diagonalization of spin chains with up to 20 
sites depending on the site spin. We have considered 201 
values of J2 in the range to 2 and 101 values of 8 in the 
range to 1 corresponding to over 20,000 grid points. 

Wc have studied the complete phase diagrams of these 
three systems using entropy contour plots and fidelity 
in the J 2 — S plane. We have been able to identify the 
quantum phase transitions from gapless to gapped phases 
using the density of the contour lines of the entropy and 
the minimum fidelity. Though the full phase diagram of 
the spin-3/2 system has not been investigated before, we 
have conjunctured the existence of some gapless regions 
and an order-disorder line by studying its phase diagram 
and comparing it with the phase diagrams of the spin-1/2 
and spin-1 systems. Our main results are that we find 
indications of gapless region near 6 = and small values 
of J2 for the spin-1/2 system, a gapless region at finite 6 
for the spin-1 system, and two gapless regions near 5 = 
and finite values of 6 respectively for the spin-3/2 system. 
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